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Abstract 

In view of the intriguing, preliminary search results for instanton-induced events at 
HERA from the HI collaboration, some important remaining theoretical issues are dis- 
cussed. Notably, the question is addressed, to which extent the HI analysis may be directly 
compared to our original predictions from instanton-perturbation theory, since certain fidu- 
cial cuts are lacking in the HI data. Various theoretical uncertainties are evaluated and their 
impact on the observed excess is discussed. An improved understanding of the experimental 
findings along with an encouraging over-all agreement with our original predictions seems 
to emerge. 



1. Instantons represent a basic non-perturbative aspect of QCD, theoretically discovered and first 
studied by Belavin et al. pQ and 't Hooft j2], about 25 years ago. As topologically non-trivial 
fluctuations of the gauge fields with a typical size of 0.3 -r- 0.5 fm [3J IH EJ, instantons play an 
important role in the transition region between a partonic and a hadronic description of strong 
interactions [B]. Yet, despite substantial theoretical evidence for the importance of instantons 
in chiral symmetry breaking and hadron spectroscopy, their direct experimental verification is 
lacking until now. 

It turns out, however, that a characteristic short distance manifestation of instantons can be 
exploited [7] for an experimental search: Instantons induce certain (hard) processes that are 
forbidden in usual perturbative QCD. These involve all (light) quark flavours democratically 
along with a violation of chirality, in accord with the general chiral anomaly relation [2] . 

Deep-inelastic scattering (DIS) at HERA offers a unique opportunity [Zj to discover these hard 
processes induced by QCD-instantons. It is of particular importance that a theoretical prediction 
of both the corresponding rate [HJ IH lTiT)] and the characteristic event signature HU [TJ] is 
possible in this hard scattering regimq^J. The instanton-induced cross section turns out to be in a 
measurable range 011011111, ^hera ~ C(30 100) pb, depending on cuts. Crucial information 
on the region of validity for this important result, based on instanton-perturbation theory, comes 
from a recent high-quality lattice simulation The main event signatures comprise a "fireball"- 
like final state with a very high number of hadrons, including K mesons and A hyperons, as well 
as a high total transverse energjU. With the help of the Monte Carlo generator QCDINS for 
QCD-instanton-induced events in deep-inelastic scattering [T2], the HI and ZEUS experiments 
at HERA are actively searching for signatures of instantons in the hadronic final state. The 
challenging experimental task is to distinguish the instanton-induced signal from the normal DIS 
final state near the edge of the available phase space. This effort seems however well worthwhile, 
since an experimental verification of such a novel, non-perturbative manifestation of QCD would 
be of basic significance. 

At the recent DIS2000 and ICHEP2000 conferences, the HI collaboration has reported preliminary 
results of a first dedicated search for instanton-induced events at HERA JT3HE1- In a phase space 
region, where a reduction of the normal DIS (nDIS) background to the percent level is achieved 
according to standard Monte Carlo models, a (statistically) significant excess of events was found 
in the HI data. While its size is at a level still comparable to the differences among standard DIS 
event generators, it is - for the discriminating observables - qualitatively similar to the expected 
instanton signal. The results presented are quite intriguing and encouraging, although far from 
being conclusive. Yet, they strongly enhance the motivation for looking more closely at some 
remaining theoretical issues related to our original predictions. 

2. Let us start off by briefly summarizing the strategy and essential results of the HI search for 
instanton-induced events jlHjllf)]. Structure and kinematical variables of the dominant instanton- 
induced process in deep-inelastic scattering are displayed in Fig. ^ Of particular importance will 
be the Bjorken- variables (Q' 2 , x') of the so-called instanton subprocess q' + g =>• I =^ X. 

The HI analysis is based on the strategy ^1] of isolating an "instanton-enriched" data sample by 

^or an exploratory calculation of the instanton contribution to the gluon-structure function, see Ref. 
2 A more extensive introduction to instanton-induced events in DIS may e.g. be found in Ref. |14j . 
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DIS variables: 

S= (e + P) 2 
Q 2 = -q 2 = -(e-e') 2 
x Bj = Q 2 / (2P ■ q) 

W 2 = {q + p ) 2 = Q 2 (1/xB ._ 1) 

x = Q 2 / (2p • q) = x Bj / z 

Variables of instanton-subprocess: 

Q>2 = - q >2 = -( q -k) 2 
x' = Q' 2 / (2g-q') 
W 2 = (q' + g) 2 =Q' 2 (l/x'-l) 

Figure 1: Structure and kinematic variables of the dominant instanton-induced process in deep- 
inelastic scattering. 

means of suitable cuts to a set of three sensitive, discriminating observables, defined in Table ^ 
(left). Three different cut-scenarios A), B) and C) with increasing instanton-separation power 
eiAnDis were defined according to 

A) the highest instanton efficiency (e{); 

B) high e\ at reasonable normal DIS background reduction (normal DIS efficiency e n Dis); 

C) highest instanton-separation power at ej ~ 10%. 

In all three cut-scenarios, significantly more events are observed than expected by standard DIS 
Monte Carlo models. It is noteworthy that with increasing instanton-separation power an increas- 
ingly large excess is seen in the data. The HI collaboration has compared the observed excess 



Observables 


^ch.b = number of charged particles 
in "instanton band" 


E% = total transverse energy 
in "instanton band" 


Sph = sphericity in the rest system of 
particles not from current jet 


1 — £out,b/-Kin,b = Et- weighted azimuthal 
event isotropy 


Q' r 2 c = reconstructed virtuality 

of quark entering J-subprocess 


Etj Ct = transverse energy of 
current jet 



Table 1: Set of discriminating instanton-sensitive observables as used by the HI collaboration ^SJ 
IT6] . Instanton-enhancement cuts are applied to the three observables on the left. All observables 
refer to the hadronic CMS {q+P = 0), except for the sphericity that is calculated in the rest system 
of the particles not belonging to the current jet. The latter is identified with the jet of highest Et. 
The "instanton band" is a prominent theoretical prediction [3 and is experimentally placed 
at fj ± 1.1 where rj is the Et- weighted mean pseudorapidity of all particles not belonging to the 
current jet. 
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in each of the six observables (Table QJ to our predictions from the QCDINS 2.0 Monte Carlo 
generator It is found that the shape and the size of the excess in the three discriminating 
observables n ch b , Sph and Q[ 2 C with instanton-enhancement cuts, as well as in the azimuthal 
event isotropy 1 — -E'out,b/-E , m,b, are in qualitative agreement with the expected instanton signal. 
However, the comparison with the two remaining observables, E% and Et-^, is more involved. 
While there appears to be also a marked excess in the respective data, both (Et^) and (Et- ]e t) 
appear to be shifted towards smaller values than predicted by QCDINS, and the widths of the 
experimental distributions are also considerably narrower. 

Figure 121 may serve as an illustration of the achievable instanton signature in the set of six 
instanton-sensitive observables (c.f. TableHJ) with cuts, integrated luminosity (/ Cdt = 15.8 pb _1 ) 
and parton densities [E] as for the preliminary HI data fSlEj- While refraining to show the real 
data before their final publication, we have instead displayed the statistical errors as calculated 
by means of the HERWIG 5.9 [TH] event generator for vanishing excess. On the one hand, the 
corresponding number of 361 ± 26 nDIS events from HERWIG after instanton-enhancement cuts 
(cut-scenario C [HHHS]) matches quite well the values quoted by the HI collaboration. On the 
other hand, the detector corrections (not accounted for here) seem to significantly reduce the 
number of instanton events by almost a factor of two. 

A central issue to be discussed in this letter is the question to which extent the HI analysis may 
be directly compared to our theoretical predictions, since certain fiducial cuts are lacking in the 
HI data: Because of difficulties with the x'-reconst ruction, there is presently no cut on x', while 
the data are compared to QCDINS with the default x'-cut implemented. Moreover, both the data 
and the QCDINS results used are lacking the theoretically required cut in Q 2 > 0(100) GeV 2 . 
We also use this opportunity to reconsider in some detail the theoretical uncertainties due to the 
extraction of the fiducial cuts in (Q' 2 ,x') from lattice data and due to the known uncertainties 
5A^jg of the QCD scale. 

3. Proceeding in reverse order for reasons of presentation, let us turn first to the uncertainties 
arising from 5A^g via the known strong dependence [5] of the instanton subprocess cross section 
(c. f. Fig. |TJ) on the QCD scale parameter, 



The dependence of the various observables on the uncertainty of is displayed in Fig. |21 The 

solid line corresponds to the default prediction of our QCDINS event generator 12J, with A™ = 
346^9 MeV, being obtained by "flavour reduction" with 3-loop accuracy from the 1998 world- 
average value PH = MeV or a s (M z ) = 0.119 ± 0.002. Note that the errors quoted 

for the total expected instanton-event numbers in Fig. |2] only refer to these A^ uncertainties. 
The optimization of the instanton-enhancement cuts (cut-scenario C of Ref. [TH| I16j) reflects in 

(5) 

the number of instanton events being maximal for the default value of A^. 

Apparently, the actual dependence of the various observables on <5A^ is much less dramatic than 
naively expected from Eq. (JTJ). The reason is, of course, associated with the fact that the fiducial 
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Figure 2: Predicted excess due to instanton-induced events (QCDINS 2.0 [121), relative to the 
normal DIS expectations (nDIS). Except for missing detector corrections, all cuts, luminosity (c. f. 
Table above), observables (c. f. TableQ) and the parton densities ^7] are as for the preliminary HI 
data ^3 ■ While refraining to show the real data before their final publication, we have instead 
displayed the statistical errors as calculated by means of the HERWIG 5.9 fXH| event generator 
for vanishing excess. Apart from modelling the achievable significance of the instanton signal via 
a good approximation of the actual (statistical) errors, the figure displays the dependence of the 

(3) 

various observables on the uncertainty of A^L The solid, dotted and dashed lines refer to the 



1998 world average value AjjJL = 219 MeV [19j, its upper and lower la limits, respectively. 



MS 



5 




p(°)(Q',...)[fm] R*/P*(x',...) 

Figure 3: Top: Illustration of the agreement of recent high-quality lattice data |3J HHI for the 
instanton-size distribution and the //-distance distribution with the predictions from instanton- 
perturbation theory P2j for p<0.35 fm and R/p> 1.05, respectively. Bottom: Display of the 
instanton-induced event distributions for HERA from QCDINS, as function of the saddle-point 
values (p*(Q', . . .), R*/p*(x', . . .)), with the default (Q' 2 , x')-cuts implemented as in Fig. |21 In 
these variables, the above cut-off values (p^L' (R/ p)mm) extracted from the lattice are apparent. 



Q'-cut in Fig. |21 is proportional to Ai^. This is a natural consequence from exploiting lattice 
results to which we shall turn next. 

4- Let us reconsider j^J ED] in more detail the extraction of the fiducial cuts in (Q' 2 ,x') from a 
recent high-quality lattice simulation of quenched QCD j^j, with special emphasis on the major 
associated uncertaintie^E 

The validity of instanton-perturbation theory, on which our predictions for HERA are based, 
requires instantons of small enough size p < p max along with a sufficiently large separation RJ p > 
(R/ p) m i n between them. Crucial information on (p max , (R/ p) m in) was obtained [H] by confronting 
the predictions of instanton-perturbation theory with the UKQCD lattice "data" (Sj for the /-size 
distribution and the //-distance distribution d4x (c.f. Fig. El (top)). 

3 A proper account of the considerable systematic uncertainties associated with different cooling/smoothing 
methods in the lattice simulations |3J|S] and 0] is very hard and clearly beyond the scope of this paper. 
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Let us recall [01 ED] that the instanton-induced cross section 0j^ RA depends strongly on these 
lattice observables via Fourier-type transformation^, 



^(<3'V) ~ / d A R e i{p+q ' ) R J dp J dp D(p) D(p) J dUe~^ n y^' Ep,U ) " «'< ~ 



e -Q'{P+ P ) | y ( 2 ) 



duj 

d^x dp 



d 4 x d 4 R 



Due to this structure, the limits (p ma x, (R/ p)min) for the validity of instanton-perturbation theory 
from lattice data may, in fact, be translated in a one-to-one manner, via a saddle-point relation, to 
minimally required cuts on the conjugate Bjorken variables (Q' 2 ,x') of the instanton-subprocess 
as follows, 

a(*\f) a(0) 

(hi 1 go_ 1 ■ 

(n f ) (n f ) ~ .(0) (0) ) 



& mm 

.(«/) j 

— I saddle-point 



MS 
„/ 



X 
mm 



(Bs.y n,) ~ (Ry 

V P* / min ~ VP/] 



According to the known n/- dependence of the solutions (p*, i?*) of the saddle-point equations [0], 
one finds that the combinations (p* A^/ /? , -R*/ p*) are approximately n/-independent functions of 
(Q'/Ajyfg, x'). In addition, Eq. Q incorporates the working hypothesis that the ^-independence 
of this combination extends also to the corresponding limits of validity of instanton-perturbation 
theory. With high-quality lattice data being only available for nf — at present, this appears to 
be a near-at-hand prescription to account for effects of light flavours in the required fiducial cuts. 
In practice, one will identify the actual (Q',x r ) cuts with the "fiducial" ones, i.e. (Q' mm , x' min ), in 
order to profit from a high event rate without loosing theoretical control. 

Before turning to a discussion of uncertainties associated with this procedure, let us briefly re- 
capitulate the main result on the fiducial region from Ref. ^0] and its implementation in the 
QCDINS event generator for HERA. First of all, Fig. El (top) illustrates the good agreement 
of the data for the instanton-size distribution and the //-distance distribution from Refs. [3 EI] 
with the predictions from instanton-perturbation theory from Refs. [TQJQ21, for 

implying J ^ ^ ms ~ ' u\ 
x' > 0.35. 

Figure El (bottom) displays the instanton-induced event distribution for HERA from QCDINS, 
as function of the saddle-point values (p^\Q', . . .), R*/p*(x', . . .)), with default fiducial cuts in 
(Q' 2 ,x') as in Fig. [21 (top). Apparently, the event distribution correctly reflects the above cut-off 
values (p^ax; (R/p) min) extracted from the lattice. 

The next issue is to extract a more quantitative estimate for the associated uncertainties. Starting 
with the p-dependence in the /-size distribution, we also take into account the error on the QCD- 
scale parameter for rif = from Ref. [201 as follows. We treat A^ as a free parameter 




MS J f-^* — —«-■■- - MS 

4 In Eq. (J2J), D, and {. . .} denote the known theoretical expressions for the /-size distribution, the //-valley 
interaction and further smooth factors, respectively. 
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Figure 4: x 2 /(# d.o.f.), characterizing the deviation of instanton-perturbation theory from the 
UKQCD lattice data CHj for the /-size distribution according to Eq. (JHJ). Left: Contour plot 
for a range of maximal values of A^p and of the renormalization-scale factor s. Darker shading 
corresponds to increasing x 2 1 (# d.o.f .). Right: Projection on the p-axis for a range of larger 



s- values. The solid line corresponds to the default A^p cut-off in QCDINS 2.0 and the dotted 
one to an estimate of the associated uncertainty. 



fluctuating around its quoted mean value, (AgQ = 238 MeV, with error 

5A^> = 19 MeV, in the 

IVlo IVlo 

following x function, 



The theoretical expression for the /-size distribution Di^ pert .th.{Pii s > A^r) from instanton-pertur- 
bation theory is taken as in Ref. ^U]. It incorporates renormalization-group invariance at the 
two-loop level, d \n(D(p, s, A^))/ d ln(s) = 0(a 2 ), with the renormalization-scale dependence 
parametrized conveniently as 

H = s/p. (6) 

The resulting xV(# d.o.f.), minimized with respect to Aj^L is displayed in Fig. 0] for a range of 

maximal values of A^p and of the renormalization-scale factor s. First of all, it is important 
to consider sufficiently large s > 1, where the incorporated two-loop renormalization-group in- 
variance should lead to a largely s-independent limit p^ x . Indeed, this is nicely reflected in the 
contour-plot of the resulting xV(# d.o.f.) versus s and A^p in Fig. HI (left). The projection on 
the A^ p-axis for 1.65 < s < 3.5 in Fig. 0] (right) indicates a fairly sharp increase of x 2 1 d.o.f.) 
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around A^p m 0.42, corresponding to our previous default cut-off p^ x ~ 0.35 fm. As a semi- 
quantitative upper limit for p^^, one may use the value, where the probability corresponding to 
X 2 /(# d.o.f .) has dropped frorrjj (9(20)% to 5%. In this way, we find 

A^pL°L = 0.42-0.435, implying ^0 = 28.2 -30.8. (7) 

MS 

As to the _R/p-dependence of the //-distance distribution in Fig. El (top, right), a naked-eye 
estimate for the uncertainty in (-R/p) mm gives the range 

= 1.0 -=- 1.05, implying x' min = 0.31 -=- 0.35. (8) 

in 

More quantitative methods do not seem worthwhile here, due to remaining theoretical ambigu- 
ities jlOj : The integrations over the /- and /-sizes in the //-distance distribution (c.f. Eq. Q) 
imply significant contributions also from larger instantons with 0.35 fm<p,p<0.6 fm, say. 

5. Let us now turn to the crucial question to which extent the HI analysis may be directly 
compared to our theoretical predictions. As mentioned earlier, possibly important fiducial cuts 
are lacking in the HI data. Firstly, we address the fact that there is presently no cut on x' in the 
data, while these are compared to our QCDINS predictions with the default x'-cut implemented. 

A glance at the lattice data for the //-distance distribution (Fig. El (top, right)) reveals that 
actually, instanton-effects seem to be very strongly suppressed, as soon as instantons and anti- 
instantons start "touching" each other, i.e. for small enough separation, R/p < 0(1). This 
rapid onset of instanton suppression corresponds to x' < 0.35, which happens to coincide with our 
default x'-cut. Hence, since our QCDINS predictions, with the x'-cut implemented, model quite 
well the actual suppression (c.f. Fig. El (top, right) and (bottom, right)), we consider the lacking 
experimental reconstruction of x' not a serious problem in principle. 

Quantitatively speaking, however, the lacking experimental x'-cut leads to a quite substantial 
ambiguity in the expected size of the instanton signal in the various observables. In this case, the 
onset of the actual instanton suppression from Fig. El (top, right), modelled by an x'-cut within 
the allowed theoretical x min -uncertainty window (JH1), gives rise to a possible increase of the size of 
the instanton signal up to a factor of three compared to our default prediction. This is displayed 
in Fig. El (left). Yet, the predicted shapes of the six distributions in Fig. El are virtually unaffected 
by this important remaining uncertainty. Upon variation of x min within the window (jHJ), we note 
that to very good approximation, the predicted instanton-induced excess in all six HI observables 
of Fig. EJvaries by the common factor displayed in Fig. El (left) for cut-scenario C. In particular, a 
slight reduction, x min = 0.325, compared to the QCDINS default x'-cut, leads to an increase of the 
instanton signal by a factor of two, representing a remarkably good description of the observed 
excess in four (n c h,b, Sph, Q' r ^ c and 1 — E out jh / E in ^) out of the six experimentally considered 
observables. While the size of the observed excess in the remaining two observables Et^ and 
Etj Ct is quite satisfactorily described as well, both (Etb) and (Etj et ) appear to be shifted towards 
smaller values than predicted by QCDINS, and the widths of the experimental distributions are 
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Figure 5: Left: Common enhancement factor for the instanton signal in all six HI observables 
(Tabled Fig- El cut-scenario C), arising upon variation of the x'-cut within the allowed uncer- 
tainty window (JHJ). Middle: Illustration of the strong x — Q 2 -correlation according to QCDINS 
for cut-scenario C. The strong concentration of events around (x ~ 0.035, Q 2 ~ 20 GeV 2 ) is 
apparent. Right: Predicted excess of charged kaons due to instanton-induced events according to 
QCDINS relative to the normal DIS expectations (nDIS). The meaning of the various lines, 
the statistical errors, cuts and luminosity are as in Fig. 



also considerably narrower. 

We shall argue next that this apparent discrepancy is actually a consequence of the experimentally 
lacking but theoretically required cut in Q 2 , 

Q 2 > Ql m « QZn = ^(lOO) GeV 2 . (9) 

As a brief reminder (HI E] , this cut assures in particular the dominance of "planar" handbag- type 
graphs in ct^era an d all final-state observables. The non-planar contributions do not share the 
simple, probabilistic interpretation of the planar ones, involve instantons with a size determined 
by 1/Q rather than 1/Q' and are both hard to calculate and hard to implement in a Monte Carlo 
generator. On account of their known power suppression in 1/Q 2 and a cross-check in the simplest 
case without final-state gluons JS], they can be safely neglected upon application of the cut 0. 
Because of these reasons, the non-planar contributions are not implemented in QCDINS. 

In order to gain some further insight as to the influence of this Q 2 -cut, consider Etj et , for example, 
for which one may easily derive the following simple expression on the parton level in the hadronic 
CMS, 

^«7^!f-^f, no) 

involving besides Q 2 , Q' 2 and x', the Bjorken-x of the 7* g =^> X subprocess in Fig. ^ The 
important issue to be exploited next is the strong correlation between Q 2 and x, apparent in 
Fig. El (middle). Already at this point, it is clear that via its x-dependence (fTUj). the Et- iet - 
distribution is strongly dependent on the underlying Q 2 -spectrum. In absence of the theoretically 

5 Clearly, the absolute value of the probability is not very meaningful, since, for the three lattice data sets used, 
the unknown but significant systematical errors have been ignored. 
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Figure 6: Dependence of the shape-normalized distributions of the six HI observables (Table ^ 
Fig. |2J) for cut-scenario C on the low-Q 2 portion of the spectrum, for which QCDINS is not 
reliable. The shaded distributions are obtained with an additional cut x > 0.15, depleting the low- 
Q 2 regime. Note the resulting substantial shift and narrowing of the -Etjet and Et^ distributions, 
with the remaining ones being largely unaffected. 



required cut in Q 2 , the ^-distribution peaks sharply around (x) ~ 0.035 according to QCDINS 
and involves a fairly small (Q 2 ) ~ 20 4- 30 GeV 2 (c. f. Fig. 03 (middle)). In fact, this dominating 
and theoretically uncontrollable low-Q 2 portion of the spectrum is responsible for the Etj e t-peak 
being located at too large a value, since 

with {j^^^GeV*} E,J3 <«M>=G(6)GeV, (11) 

in good agreement with the direct (but unreliable) result from QCDINS in Fig. |21 not involving 
a cut in Q 2 . 

With this insight, let us perform the following exercise. We deplete the low-Q 2 portion of the 
spectrum by means of a cut, x > 0.15, in the (internal) Bjorken-x variable. The QCDINS results 
should now be about reliable since from Fig. (middle) we have effected a shift to (Q 2 ) ~ 95 GeV 2 . 
The striking result of this exercise is shown in Fig. |H1 for the six shape-normalized distributions 
of the HI observables. First of all, we observe that the shapes of the four distributions of n ch b , 
Sph, Q' T 2 C and 1 — E out h /E irijb , for which there was good agreement with the observed excess 
before applying the x-cut, remain essentially unaffected. Their shapes are thus rather insensitive 
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to the Q 2 -spectrum and should be reliably predicted by QCDINS, even without application of 
the required cut in Q 2 . On the contrary, the shape- normalized distributions of Et )et and E% 
appear both shifted significantly towards smaller values and are now much narrower. Notably, the 
shape-normalized E'tjet-distribution is now in virtually perfect agreement with the corresponding 
experimental quantity. Also the .E^b-distribution matches well the experimental trend, after this 
depletion of the low-Q 2 region. 

On the one hand, these arguments may well serve as indication that upon proper implementation 
of the required fiducial cuts the good description of the data will persist. On the other hand, 
they should be taken as a strong encouragement to try and implement the lacking cuts in Q 2 and 
x' into the data. The essential, remaining problematics is that searches for an instanton-induced 
excess in DIS have to rely on background estimates from nDIS Monte Carlo generators, in a 
region of phase space where their accuracy is not too well known. 

We close with a reminder of our predictions for the rate of charged kaons that can and should 
be investigated by means of dE/dx in the instanton-enhanced data sample (c. f. Fig. El (right) 
for cut-scenario C). An excess of kaons, as direct consequence of the flavour democracy of the 
instanton-induced interaction, represents a crucial independent signature. Respective data will 
hopefully be included in the final HI analysis. 
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